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I. INTRODUCTION 

During the last several years, there has been a revival 
of interesti*2i2iii5iSiL2i2ii2iiiii2i in unconventional types of 
spin ordering in systems with higher spins S > 1. This 
interest is to large extent motivated by experiments on 
Bose-Einstein condensates of cold atoms with internal 
spin states ! 13 ' 14 Particularly, ordering of the quadrupole 
degrees of freedom corresponds to the so-called "spin- 
nematic" type of spin order , 15 ' 16 ! 17 which is difficult to 
obtain in conventional magnetic materials since its ex- 
istence requires the presence of strong non-Heisenberg 
(biquadratic or multispin) exchange terms, or the pres- 
ence of strong frustration mixing ferro- and antiferro- 
magnetic couplings ! 18 ' 19 ' 20 In strongly frustrated sys- 
tems with ferromagnetic couplings, higher multipolar or- 
ders may win over the nematic one, becoming dominant 
correlations ! 19 ' 20 ' 21 For cold spinful bosons in optical lat- 
tices, strong non-Heisenberg exchange appears in the ef- 
fective spin model ) 22 ' 23 favoring spin-nematic order. In 
higher-spin systems, higher symmetries may naturally 
arise. SU(N) generalizations of Heisenberg spin systems 
in one and two spatial dimensions have been extensively 
studied^ 25 ' 26 ! 27 ' 28 ' 29 Several recent studies^^ 2 . ex- 
plore exotic pairing possibilities opened by the existence 
of higher SU(N) symmetries with N > 2 in fermionic 
systems. 

In the present paper, we will study what happens to 
an SU(N) antiferromagnet (AF) if the high symmetry 
gets explicitly broken by a weak perturbation. It will be 
shown that, similarly to N = 2, the physically important 
case N = 3 is in many respects special, and breaking the 
SU(3) symmetry leads to rich behavior which might be 
realizable in cold atom setups. We will see that perturb- 
ing the SU (3) symmetry has a drastic effect on the topol- 
ogy, which is reflected in physical properties due to the 
role of the Berry phases. Our starting point will be the 
5 = 1 model on an anisotropic square lattice described 
by the Hamiltonian 

ii 

h n ,n> = cos9(S n ■ S n >) + sm8(S n ■ S n >) 2 , (1) 

where S n is a spin-1 operator at the lattice site n, and A 
is the parameter controlling anisotropy of the lattice, < 



A < 1. This model appears, particularly, in the physics 
of ultracold alkali atoms with hyperfine S = 1 spins (e.g., 
23 Na) in optical lattices at odd filling.^2, The parameter 
6 can be varied by tuning the ratio 02/00 of scattering 
lengths in S = 2 and S = channels using the Feshbach 
resonance, as well as by creating a gradient in the optical 
lattice potential. 33 Similar models have been proposed 7,8 
as a possible explanation for the unconventional spin 
state discovered recently^ in the quasi-2d S = 1 mag- 
net NiGa2S4, and have been also discussed 9 "^ in context 
of the deconfined quantum criticality conjecture 3 ^. In 
one dimension (d = 1), this model has been extensively 
studied and a number of analytica l 35 ' 36 '^ 7 .' 3 ^' 3 ^ -'^- and 
numericalMi 5 -^^ 2 - results are available. In two dimen- 
sions, it was recently studied numerically by means of 
Quantum Monte Carlo technique 9 and analytically with 
the help of a field-theoretical approach] 10 ' 43 

Using the standard representation of the S = 1 opera- 
tor S% = —it a p 1 ^ n pt nrf through three bosonic operators 
t a , a = 1, ... 3 satisfying the hardcore constraint 

tU a =n c = l, (2) 

one can cast the local Hamiltonian in the form 

h ■ ■ — — Tt^ ft t t —Ji^i^i t 

J = — cos 8, J = cos — sin 9. (3) 

Since the model is formulated in terms of local bilinears 
of bosonic operators t n , it obviously has the local U(l) 
symmetry for any values of the model parameters. We 
will be interested in the interval — 3ir/4 < 9 < 0. It is 
convenient to generalize the Hamiltonian by letting 
the boson flavor index run from 1 to TV and allowing the 
parameter n c in ^ to be an arbitrary integer number. 
In case of the related models for cold atoms in optical 
lattices, n c has the meaning of the number of atoms per 
lattice site^ and in what follows we will assume n c to 
be odd. For n c = 1, N = 3 corresponds to |T]), N = 2 
describes the spin-i XXZ model with J x = — 2(J + J), 

J z = 2(J — J), and N — 4 can be realized^ as a "bi- 
layer" spin-i model with four-spin interaction between 
the layers: 

h n , n ' = (2cos0-sin0)[(s n -s n /) + (T n -T n ')] 

+ 4 sin 9(s n ■ 8 n /)(r n • T n >), (4) 
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which is essentially the Kugel-Khomskii spin-orbital 
model 4 ^ with spin and orbital degrees of freedom de- 
scribed by s n and r„ spin-^ operators, respectively; a 
large number of results are available for this model in 
one dimension^^^^^^^ 4 .^ 

The point 9 = —tt/2 (J = 0) is remarkable since it 
has an enhanced symmetry. The Hamiltonian ([1]) is al- 
ways SU(2) (generally, O(N)) invariant, but at 6 = —tt/2 
the symmetry group is enlarged to SU(3) (respectively 
SU(N)): since the lattice is bipartite, a transformation 
t n i ► Ut n on sites n belonging to sublattice A leaves 
the Hamiltonian invariant if it is accompanied by a con- 
jugate transformation t n t— > U*t n for n G B, with a 
unitary matrix U . 

Our strategy will be to construct an effective field- 
theoretical description of the problem, using 9 = —tt/2 
(J = 0) as a starting point, and to treat the term pro- 
portional to J as a perturbation. We will also see that 
a rich behavior is generated by adding another pertur- 
bation, namely the easy-axis single-ion anisotropy to the 
5 = 1 Hamiltonian |T]), 



h^h-dY,(s z J 2 , d>o. 



(5) 



For the case of a general N this amounts to including 
the term of the form Dj^ n t n N^n N> which breaks the 
symmetry down from SU(N) to SU(N — 1). In cold 
atom systems, such terms appear naturally in presence 
of external magnetic field due to the quadratic Zeeman 
effect 

The structure of the paper is as follows: in Sect. ITT1 the 
effective continuum theory in the vicinity of the SU (N)- 
symmetric point is derived, Sect. ITTT1 considers the influ- 
ence of the 9 perturbation breaking the symmetry down 
to O(N), Sect. IIVI studies the effects of the anisotropy 
flSJ, and, finally, Sect. IVl contains a brief summary. 



II. EFFECTIVE FIELD THEORY 
IN THE VICINITY OF THE SU(N) POINT 

To construct the continuum field description, consider 
a path-integral representation of the problem, effectively 
replacing the bosonic operators t n>a with complex fields 
on the lattice satisfying the constraint in (J2J) . To pass to 
the continuum properly, one should notice that local spin- 
quadrupolar correlations are of the ferromagnetic type 
for J + J > 0, while the spin-dipolar correlations are 
antifcrromagnetic provided J > Jj££ this can be also seen 
from the numerical results (see Fig. 8 of Ref. |42) . We 
will be interested mainly in the region of 9 < 0, where 
the first of those inequalities is always satisfied, but the 
second one breaks for 9 < 9q w — 0.657T. The theory 
derived here will be valid for 9 > 9q, and the proper 
effective theory for 9 < 8q can be found in Ref. l58l^ 

The AF character of local spin correlations suggests 



the following ansatz for the bosonic lattice fields t n : 

tn = K + iVnVn) + {Vn<fn + %Cn), (6) 

where rj n is an oscillating factor taking value ±1 for n 
belonging to A and B sublattices, respectively, and u, 
v, ip, and C are assumed to be smooth functions of the 
site coordinate n. Defining z n = (u n + iv n )/^/n^ and 
ip n = (<p n + i£ n ) I yjn~c.i one can rewrite the above ansatz 
in a simpler form: 



where the constraints 



z n + i> n , neA 

Z*n-i>*n, neB 



\ip\ 2 = 1, ip-z* +z-tp* = 



(7) 



(8) 



are implied. One can expect that the magnitude of t/>, 
which corresponds to ferromagnetic fluctuations, will be 
much smaller than that of z. Using the ansatz ([7|), pass- 
ing to the continuum, retaining only up to quadratic 
terms in ip and neglecting its derivatives, one readily ob- 
tains the Euclidean action A = Ao + Amt + Ab, where 
Aq corresponds to J = 0: 

Aq = V\n 2 c I ' dr J d 2 x^{i>* ■ d T z - ip ■ d T z*) 



+ 4/(1 + A) 
+ J[\d, 



■ d k z\ 2 ] 



(9) 



+ //!(</> ■ + r ■ z) + MN 2 + M 2 - 1)}, 



and Ab is the topological Berry phase contribution 

A B = in c ^2 r} n arg (z* (r) • z n (r + dr)) (10) 



which is known to play a crucial role in the physics of 
the systemi 24 i 25 i 60 It is important to realiz o 10 ! 22 that the 
naive continuum limit of (|10p . 



A B = n c J dr^2 VnZ n ■ d T z n , 



(11) 



can only capture the contributions from smooth field con- 
figurations and in case of dominant nematic correlations, 
when z becomes a real vector defined up to a sign, misses 
the additional phase stemming from disclinations. 

The term Aint is determined by the "perturbation" J, 

Aim = Jy/Xnl Jdr J d 2 x{-\z 2 \ 2 + \z-d k z\ 2 

+ [z 2 r 2 + \z\d k z*) 2 + c.c]). (12) 

Here the index k runs over two spatial coordinates, the 
factor VA in ([5]), ([12")) comes from rescaling one of those 



3 



coordinates to compensate for the anisotropy of interac- 
tions, and fj,i 2 are the Lagrange multipliers ensuring the 
constraints. 

For J — 0, one can easily integrate out tp and /ii fields; 
it turns out that /ii = — n~ 1 z* ■ d T z, which yields 

V> = -[4Jn e {l + \)]- 1 {a r z-z(z* -d T z)}, 
= [Un c {\ + X)]- 1 {d T z* + z*(z* -d T z)}. (13) 

Substituting this back into ([9]), one obtains the effective 
action for z field only, where we can now approximately 
assume \z\ 2 — 1. Rescaling the imaginary time axis r i— > 
t I (2n c J\J\ + A), one arrives at the effective action 

A ° = Ygj d d+1 x{\d»z\ 2 -\z*-d„z\*}, g = ^ + A "' , 

^(14) 

where d = 2 is the spatial dimension, and the index /i 
runs over all d+1 space-time coordinates. Had we started 
with a single 5=1 chain instead of the square lattice, we 
would have obtained the action of the same form (fl4|) , but 
with d — 1 and g = l/n c . This is nothing but the action 
of the CPN" 1 model^i^SiSlSa originally proposed as 
an effective theory for SU(N) antiferromagnets by Read 
and Sachde v 24 ' 25 . This action has a local U(l) gauge 
symmetry z i— ► e^^z and can be rewritten in the form 

A ° = ^J d d+1 x\(d„ - iA^zl 2 , (15) 

where = i(z ■ d^z*) is the U(l) gauge field. 

The CP 1 ^^ 1 model without the topological phase term 
is always gapped in d = 1 and displays an ordering tran- 
sition in d = 2 at a certain critical value of the coupling 
constantj 63 ' 64 In the disordered phase the z field acquires 
a finite mass, and a kinetic term for the gauge field is dy- 
namically generated,— 

A^A+^J d d+1 xF 2 ^ (16) 

where F^ v = <9 M A„ — <9„A M , and the coupling constant 
e 2 cxA 3 - d . 

The Berry phase term (|T0|) is crucial for the physics 
of the disordered phase ; 24 ' 25 ! 60 particularly, it leads to 
spontaneous dimerization in d = 1 for odd n c (except 
for N — 2 which is special: in that case the system re- 
mains gapless and translationally invariant in a wide g 
rang o 66 ' 67 ' 68 ' 69 ), and in two dimensions the disordered 
phase gets spontaneously dimerized in different patterns 
depending on the value of (n c mod 4). We will come 
back to the role of the Berry term later and look into the 
rest of the action first. 



III. EFFECT OF THE SU(N) ^ O(N) 
PERTURBATION 

The perturbing action (|12[) explicitly breaks the global 
SU(N) symmetry down to O(N), but preserves the U(l) 



gauge symmetry. Consequently, nonzero J can produce 
only gauge-invariant perturbation terms of the form 

|* 2 | 2 , \z-D^z\\ 

where D M = <9 M — iA^ and the ellipsis stands for terms 
with higher derivatives. It is easy to see that the first 
term above is relevant for d < 3, while the second one 
is irrelevant for d > 1 (for d = 1 it is marginal). It 
thus makes sense to consider only the effect of the most 
relevant term, which brings us to the perturbed action 
Aj = Ao + Ant , with 

Ant = -^/d d+1 *NT, 7- J- (17) 

It is easy to generalize the standard large- N mean-field 
analysis 63 of the CP N_1 model to include the effect of the 
SU(N)-breaking perturbation 7. We consider the action 

Amf = Y g J d d+1 x{\d,z\ 2 - \z* ■ d^z\ 2 

- 7 |* 2 | 2 +a(N 2 -l)}, (18) 

where a is the Lagrange multiplier responsible for the 
constraint \z\ 2 — 1, and expand it around a stationary 
saddle-point solution z = Zq, a = <jq. This expansion 
has to be performed differently depending on whether 
the perturbation is of the "nematic" (7 > 0) or "antifer- 
romagnetic" (7 < 0) type. 

A. "Nematic" side (7 > 0) 

In this case the saddle point can be chosen in the 
form Zq = (no, 0, . . . , 0), which in our original N = 3 
model corresponds to the spin-nematic (quadrupolar) or- 
der. Fluctuations around the mean-field solution, z = 
z + u + iv can be described by two real TV-component 
vectors u, v. Due to the constraint \z\ 2 = 1 one can set 
Mi w 0, and the gauge-fixing condition (e.g., setting no to 
be real) yields »i » 0, After integration over quadratic 
fluctuations, the saddle point equations are obtained as 

n 2 +g(jv-l)^{^-+ 1 ^, 2 } = L 

h 

i no = 0, (19) 

where the sum is over (d + l)-dimensional reciprocal 
space. In one spatial dimension, d = 1, the model is 
disordered (no = 0) for any value of the coupling con- 
stant g, and the field z is always massive, with ao — A 2 
having the meaning of the squared spectral gap, 

A 2 cA 2 exp{--^ Ty }-2 7 , 7 «A 2 , (20) 

where A is the lattice (UV) cutoff and it is assumed that 
7 <C A 2 and both 7 and A are small compared to the 
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FIG. 1: A sketch of the phase diagram of the model ((3| on 
an anisotropic square lattice in the vicinity of the SU(N)- 
symmetric point — —n/2. For N < N c ~ 5 the phase 
boundary has a square-root- type cusp at 6 — > — tt/2 as sug- 
gested by Eqs. l|2"2|). (f2"o) . For N > N c the dimerized phase 
has a finite extent at A = 1. The phase denoted as "Haldane" 
is for N — 3 indeed the Haldane phase whose boundary lies 
at 6 — — 7r/4, and in case of = 4 it is the staggered dimer 
phase^^ with the boundary at = 0. For N = 3 only, the 
degeneracy of the dimerized phase is twofold for 8 > —n/2 
and fourfold for 9 < -tt/2, see Sect. lHICl 



cutoff. In the opposite case if 7 ^> A 2 one obtains 



A A2 



{-m^}> - >>A2 (21) 



In two dimensions, d — 2, there is a finite second-order 
transition point g = g c , given by 



N 



rV7), 



(22) 



such that for g < g c the O(N) symmetry is sponta- 
neously broken and the ground state is ordered, n 2 , = 
1 — g/g c , while for g > g c one has a disordered phase 
with n = and <7rj = A 2 , where the gap A behaves as 
A^ jtiidc 1 -^ 1 ) at 5 - 5c . 

The transition at g — g c corresponds in our origi- 
nal model ([I]) to a transition at some critical value of 
anisotropic coupling A = A c , so that one has the spin- 
nematic ordered phase at A > A c and the quantum dis- 
ordered phase at A < A c , and the "disordered" phase 
actually corresponds to a dimerized state arising due to 
the Berry phase term^ The critical value A c can be esti- 
mated using known large- N result* 2 ^ for the critical point 
at 7 = and the isotropic square lattice (A = 1), n™' w 
0.19AT, which yields A" 1 (7 = 0) w 55A(n c /N) 2 - 1. Since 
only < A c < 1 makes sense, the latter estimate suggests 
that in absence of the perturbation 7 the system does not 
order for any value of A for N > N c ~ 5. Extrapolating 
to N = 3, one obtains that at 7 = (i.e., 9 = — tt/2) the 
critical coupling is A c w 0.19, while the Quantum Monte 
Carlo calculations^ done for bilinear-biquadratic 5 = 1 
model yield A c ~ 0.13. 



"Antiferromagnetic" side (7 < 0) 



For 7 < the interaction favors the minimal absolute 
value of z 2 , so the saddle-point solution can be chosen in 
the form z$ — 2 _1 / 2 (no, ino, . . . , 0), which in the original 
N = 3 model corresponds to a finite AF order parameter 
£ = —i{zq x Zq). For the fluctuations it, v one can set 
U2 ~ 0, v\ fa to fix the gauge, and the constraint 



\z\ = 1 yields u\ 
take the form 



'-'2 



0. The mean field equations 



N — 2 
ao + k 2 



1 



cr +4|7| + k 2 



<7 n = 0. 



(23) 



For d — 1 there is again only a disordered phase with 
and ctq = A 2 , where 



n 



A2exp {-^T)} 



4|7l 

g(N-l)) N-l 



(24) 



under the assumption I7I -C A 2 , and in case A 2 ^> I7I 
the gap is given by 



\4h J 



exp ■ 



(25) 



*W 1 9(N-2)i 

For d = 2 the system orders at 5 below the critical 
value g c given by 



9c 1 



A(N-l) 



(26) 



and for g > g c one has a disordered phase with a finite 
gap A = y / cto which grows linearly in the vicinity of the 
transition, A ~ ^(g,: 1 - g^ 1 ). 

The corresponding phase diagrams are sketched in Fig. 
[TJ On the AF side the effect of perturbation 7 is weaker 
by a factor of ~ 1/N compared to the "nematic" case 7 > 
0: for d — 2 this is translated into different amplitudes 
of the square-root cusp in the dependence of the critical 
coupling A c on 7 for 7 > (8 < —n/2) and 7 < (9 > 
—tt/2). For d = 1 this effect should be seen in different 
slopes of the gap A(7) for positive and negative 7; this 
is in line with the results from exact diagonalization of 
small finite chains^ as well as with the recent density 
matrix renormalization group calculations for the model 
(P) on a ladder £ 



C. Influence of the perturbation on the Berry term 

Up to now we have considered only the effect of the 
S'L/(A r )-breaking perturbation 7 on the action without 
the Berry term. Apart from favoring nematic or antifer- 
romagnetic order, the effect consists in a mere shift of 
the transition point in two dimensions, and a change of 
the gap in d = 1 case. However, there is another im- 
portant effect of the perturbation 7: as we will see, it 
drastically affects the Berry term, which has important 
consequences for the physics of the disordered phase. 
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1. d=l 

The role of the Berry phase term Ab at 7 = has been 
studied in detail^^^ In the onc-dimensional case one 
obtains 



A 



(d=i) 



iQq, 



1 

2^ 



dxdrF XT , (27) 



where the integer number q has the meaning of the net 
topological charge (skyrmion number), and 

= (irn c mod 2n) 

is the so-called topological angle. Explicitly expressed 
through z, the topological charge reads 



1 



^ J d 2 xe t _ tu (d t _ t z* ■ d v z) 



(28) 



For even n c the Berry phase has no effect, while for odd 
n c it leads to the twofold degenerate ground state with a 
finite "static electric field" (topological charge density) 



(iF XT ) = ± 



TV 



(29) 



One can easily show that the topological charge density 
is directly proportional to the dimerization order param- 
eter, in essentially the same way as it has been done 7 -^ 
for the 0(3) nonlinear sigma model. Indeed, the dimer- 
ization operator at 7 = can be defined as 



Ofr = Vn(S n ■ S n+1 ) 2 



(30) 



and after passing to the continuum its leading non- 
oscillating part will take the form 



h-> 2(V ■ d x z* + ip* ■ d x z) 

1 — ^ [2Jn c {\ + \)]~ 1 (d T z* ■ d x z - d T z ■ d x z*) 

= [2Jn c {l + X)]- 1 {iF XT ) (31) 

Thus, at the SLf(TV)-symmetric point 7 = the ground 
state for odd n c and TV > 3 is spontaneously dimerizedi^ 
The case TV = 2, however, is an exception: for TV = 2 the 
model is equivalent to the 0(3) nonlinear sigma model 
with the topological angle 9 = ir, which is gapless and 
nondimcrizcd 6 — i 67 ' 68 in a wide range of the coupling g. 

Let us first illustrate the effect of the SU(N) 1— > 
O(N) perturbation 7 on the Berry phase by a simple 
observation^ 7 - valid for TV = 3. Finite 7 < favors 
field configurations of the antiferromagnetic type, namely 
z = (ei + 62), with ei,2 being two orthonormal vec- 
tors and n(9, tp) — e± x e-i having the meaning of the 
unit Neel vector characterized by two spherical angles 9 
and ip. It is a straightforward exercise to check that 



2?r 



d 2 x sin Oefiu (dfj,6) (d v <p) 



= = J tfxe^n ■ (d^n x d v n) = 2Q, (32) 



where the topological charge Q is the winding number 
of the S 2 1— > S 2 mapping characterizing the space-time 
distribution of the unit vector n(9, <p). This shows that 
negative 7 favors 2-field configurations with even charge 
q and suggests that configurations with odd q become 
suppressed. This is physically important, because if odd- 
q configurations are prohibited, the Berry term obviously 
becomes ineffective, irrespectively of whether n c is even 
or odd. The above argument cannot be applied for TV > 3 
because the second homotopy group of 0(N > 3) sigma 
models is trivial so they possess no 7r 2 topological charge. 
It is also not possible to extend this argument to 7 > 0, 
because in this case "nematic" configurations with z be- 
ing a real (up to an arbitrary overall phase) unit vector 
are favored and for such configurations the CP 1 ^^ 1 topo- 



logical charge l|28|) identically vanishes. 

To understand what happens in case of general TV and 
7, consider the general one-skyrmion (q = 1) solution of 
the 1 + 1-dimensional CP N_1 model which has the form 



c a {Z — a a ) 



1/2 : 



(33) 



where Z — xq + ix\ is the complex coordinate, the com- 
plex numbers a a have the meaning of coordinates of the 
TV skyrmion constituents (sometimes called "zindons" 
from a Persian word meaning prison 7 ^), and another set 
of complex numbers c a may be viewed as amplitudes as- 
sociated with each zindon. Normalizing the amplitudes 
c a as ^ Q |c Q | 2 = 1, putting the origin into the "center of 



mass" (which amounts to demanding J2 a l c ° 



0), 



and defining the average "size" R of the skyrmion as the 
dispersion of the zindon positions, 



i? 2 = ^|c a | 2 |a Q | 



(34) 



one can recast the general q = 1 solution (f3"3")) in a more 
elegant forrn^ 



UR + VZ 



(\Z\* + R 2 



1/2' 



(35) 



where U, V are two orthonormal complex TV-component 
vectors, 



U* ■ U = V* ■ V = 1. U* ■ V = 0. 



(36) 



For 7 = 0, i.e. in the unperturbed CP 1 ^^ 1 model, the 
action of such skyrmion solution does not depend on its 
parameters. For a finite 7, however, one gets an ad- 
ditional contribution to the action from the \z 2 \ 2 term 

Let us calculate this correction to the first order in 7. 
Consider first the "antiferromagnetic" case 7 < 0. To 
minimize the action cost, we must reduce as much as 
possible the deviations of z 2 from 0. Requiring U 2 = 
ensures that z 2 — > at \Z\ — > 00, killing the next leading 
term in Z fixes U • V = 0, and, finally, if we were able to 
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satisfy additionally V 2 = 0, then the condition z 2 = 
would be identically fullfilled. Those three constraints 
can be satisfied together with ([36| only if the four real 
iV-component vectors Re(f/), Re(V), Im(C7), Im(V) are 
mutually orthogonal, which is readily achieved for N > 4 
but is obviously impossible for N — 3. Thus, for N > 4 
the q = 1 skyrmion \35\) remains an exact solution even 
for finite 7 < 0. In other words, when 7 < is switched 
on, the "zindons" constituting a skyrmion are able to 
adjust themselves for N > 4 in such a way that the 
skyrmion continues to provide the minimum of action. 
This is in fact amusing because formally for 7 < the 
model has only the O(N) symmetry, and one would ex- 
pect that skyrmions do not exist for N > 3. 

For TV = 3, 7 < the minimum contribution of the per- 
turbation to the action of the skyrmion (|35[) is achieved 
if V is real and the three vectors V , Re(Z7), and lm(U) 
are mutually orthogonal. The excess action due to finite 
7 is then given by 

AA, <0 = -^fd 2 x \z 2 \ 2 = -^, (37) 

and it grows as a square of the skyrmion size which means 
at 7 < the field configurations with q — 1 are prone to 
collapse and only exist as metastable "excitations". At 
the same time, one can easily adjust the parameters of 
a general q = 2 skyrmion solutions of the 7 = model 
so that z 2 = is identically satisfied (see Appendix |A"| . 
This effect can be interpreted as "topological pairing" 
of q = 1 skyrmions The total topological charge den- 
sity F XT can be separated into two parts, i^ir -1 ^ and 
Fxt~ 2 ^ which correspond to the contributions from un- 
bound q = 1 skyrmions and their bound pairs, respec- 
tively. Only .Fjf contributes to the nontrivial part 
of the Berry phase (|2T|) . while the full F XT enters the 
gauge field kinetic energy term (|16p . The dimerization 
order parameter (O dlm ) cx (iF XT ) will be proportional 
to the fraction p of the q = 1 skyrmions and so will be 
diminishing with increasing |7|. This "topological" sup- 
pression of dimerization at 7 < exists only for N = 3 
and is absent for N > 4, which implies that within our 
description the dimerized phase for N > 4 should extend 
to the entire region 7 < (i.e., up to 9 — ir/4 which 
in our notation is another 5?7(A r )-symmetric point of 
the model, corresponding to the transition into a criti- 
cal phased); however, it is clear that our description will 
eventually break down as the "perturbation" I7I becomes 
large. In fact, according to exact results (see Ref<22) for 
one-dimensional SO(N) generalizations of the bilinear- 
biquadratic model ([T]), even for N > 4 there is still a 
phase transition on the way from the AF SU(N) point to 
the critical SU (N) one. The chunk of different phase ly- 
ing between the dimerized phase and the critical SU (N) 
point gets squeezed with increasing N, and the transition 
point for N > 4 lies in the region of 7 > 1, way beyond 
the range of applicability of the present approach. 

Consider now the perturbation of the opposite sign 
7 > 0, which favors nematic-like field configurations with 



z = e la <p, where a is an arbitrary phase and tp is a real 
unit vector. For such z the topological charge (|2"8|) is 
identically zero, which indicates that skyrmions with any 
charge are suppressed by the perturbation. In a differ- 
ent way one can see that by calculating the 7-dependent 
correction to the action. For a q = 1 skyrmion (|35|) min- 
imizing the deviation of \z 2 \ from 1 leads to the require- 
ment that U and V are real, and the resulting correction 
diverges logarithmically with the system size L, 

. . 2irjR 2 L 

AA 7>0 ^ — —In-. (38) 
9 K 

In the disordered phase, one expects that the system 
size L above will be replaced by the correlation length 
£. In contrast to the AF-like case 7 < 0, this suppres- 
sion persists for any number of the field components N. 
A similar calculation for q = 2 yields A„4*>q cx jR 2 , 
so the even-charged skyrmions are suppressed as well, 
though weaker than the odd-charged ones. Thus, with 
increasing 7 the contribution from smooth field configu- 
rations (skyrmions) to the Berry phase dies out, but at 
the same time the contribution from discontinuous con- 
figurations (disclinations, or Z2 vortices) remains unaf- 
fected and gradually becomes the leading one. Indeed, 
a configuration with a real vector z abruptly changing 
sign across some bond along a path running in the time 
direction contributes the Berry phase equal to n for every 
such bond,— >21 which is not captured by the continuum- 
limit expression (|27p but is readily seen from the general 
formula (JTUJ) . In the disordered phase the fluctuations of 
z are gapped and can be integrated out, leaving one only 
with Ising-like degrees of freedom marking bonds where a 
discontinuous change z — > —z occurred.— The resulting 
so-called odd Z2 gauge theory^ is known to be always 
dimerized in one dimension^ which, according to Grover 
and Scnthilfi^ explains why the dimerized phase extends 
all the way up to 9 = — Sir/ 4 (their arguments can be lit- 
erally transferred to the effective theory of Ref. |H which 
is suited for describing the region — 37r/4 < 9 < ~ 0.657T 
with ferro-type local correlations). 

2. d = 2 

In two dimensions the Berry phase is determined by in- 
stanton processes ( "monopoles" ) changing the skyrmion 
topological quantum number q and is given b y 25 ' 60 

•Afl = -2 £ X)C( r ')3<. (39) 

Ti 

where the sum is over the locations r j of monopoles hav- 
ing the charge qi (i.e., the skyrmion number gets changed 
by qi), and the factor C( r i) takes values 0, 1, 2, 3 for 
belonging to the four dual sublattices W, X, Y, Z respec- 
tively (see Fig. 7 of Ref. H). At the S'J7(A r )-symmetric 
point 7 = for n c 7^ mod 4 the Berry term leads to 
the ground state with nonzero instanton density, thus to 
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finite electric fields and to spontaneous breaking of trans- 
lation symmetry^ the dimerized ground state is twofold 
degenerate for n c = 2 mod 4 and fourfold degenerate for 
n c = (1 or 3) mod 4. 

When the S'C/(iV)-breaking perturbation 7 is switched 
on, the monopoles are transformed in a similar way as in 
d = 1 case for skyrmions: at 7 < monopoles with odd q 
are strongly suppressed for N — 3 and remain unaffected 
for N > 4. Suppression of odd-charged monopoles for 
N = 3 and 7 < can be understood by invoking the 
same type of argument as that we have used in the one- 
dimensional case. The monopole charge q — § j a dS a 
can be defined as the quantized flux of the "skyrmion 
current" 



3c 



1 

2^ 



dAu 

' dx Ll 



2lT 



(40) 



through a closed surface surrounding the monopole. For 
antiferromagnetic-type configurations favored at 7 < 
one again can write z = ~^{ e i + e ■^l)^ where e% t 2 are two 
orthonormal vectors, and define the corresponding 0(3) 
unit vector field t = e\ x e 2 . Then it is easy to obtain 

j a = -^-£ af 2iy {d^e 2 ■ d u ex) 
Ztt 

= 77-£a M <4ei • (e 2 x d^)] [e 2 • (e 1 x d v ei)\ 
Ztt 

= ^-£ a ^i-{d„£xd v £) = 2J a , (41) 

where J M is the corresponding skyrmion current of the 
0(3) nonlinear sigma model whose flux through a closed 
surface should be an integer number. Again, this argu- 
ment only works for N = 3. For N > 4 the q = 1 
monopole solution 76 



z = U cos(0/2)e^ + V sin(6>/2), 



(42) 



where 9 and (p are the angular spherical coordinates in the 
(2+l)-dimensional space, and the monopole is assumed 
to be placed at the origin, can be easily adjusted to yield 
z 2 = identically for N > 4, and for N = 3 the ex- 
cess action due to the perturbation 7 of such a monopole 
diverges as the spacetime volume (note that this contri- 
bution arises due to deviation of z 2 from and thus is 
not destroyed by the vanishing spin stiffness in the disor- 
dered phase as the contribution from the main Aq part 
of the action doesZ 6 -). Even-charged monopoles can be 
shown to survive for a finite 7 < as exact solutions 
(see Appendix [A"| . So, we come to the conclusion that at 
7 < the odd- charged monopoles get confined into pairs 
for N = 3, but are insensitive to the perturbation for 
N > 4. The consequence for N ~ 3 is that the contri- 
bution of odd-charged monopoles is switched off for any 
finite 7 < 0, which effectively amounts to doubling n c in 
([55)1 : for the bilinear-biquadratic S = 1 model ([TJ) that 
means that the dimerized state is doubly degenerate at 
7 < and becomes fourfold degenerate only at 7 = 0. 



On the nematic side (7 > 0) the effective theory has 
been constructed by Grover and Senthil; 10 they have 
shown that the problem can be mapped to an XY model 
with a fourfold anisotropy term. The dimerized ground 
state is respectively predicted to be fourfold degenerate 
in that case. One is thus led to conclude that 9 = —it/2 
for N — 3 should be the first order transition line. 



IV. EFFECT OF THE SU(N) i-> SU(N - 1) 
PERTURBATION 



Consider now a different way to perturb the SU(N) 
symmetry, namely let us introduce a finite mass for one 
of the components of the z field, 



2.9 



d d+1 x\z N \ 2 , 



(43) 



which for the S = 1 model {T]) is equivalent to includ- 
ing the easy-axis single-ion anisotropy term ([5]) with 
D = 2gm\. For cold atoms in optical lattices, such 
terms appear naturally in presence of external mag- 
netic field due to the quadratic Zeeman couplingi 23 i 56 
This perturbation breaks the SU(N) symmetry of the 
model down to SU(N— 1) and produces a CP N ~ 2 model 
with the topological angle 9 — tt as the effective the- 
ory. Actually, the operators N* z = J2 n 7 lntl l _ N t n N and 

Nq z = J2 n tli N^n n commute with the Hamiltonian Q 
if cos 6* = and sin^ = 0, respectively. So, at the "ferro- 
SU(3)" point 8 — —bir/A the single-ion anisotropy D 
acts simply as an "external field" coupling to a conserved 
quantity, but at the "AF-5J7(3)" point 9 = -tt/2 the sit- 
uation is different. 

In one dimension (d = 1), if the mass m-o is large com- 
pared to the gap A ~ Aexp{— n/g(N — 1)}, one can 
integrate out just the single most massive iV-th compo- 
nent and obtain in that way a correspondence between 
the bare coupling constant g^-i of the effective CP N ~ 2 
model and the bare coupling constant g = gjy of the 
original CP 1 ^^ 1 model: 



9N-1 



1 



£M1 



(1 + ^) 

V mi / 



(44) 



Now, the case N = 3 is again exceptional because the 
CP 1 model with the topological angle 8 = tt in d = 1 
is gapless in an extended range of coupling. At infinite 
coupling 92 = 00 the parity is broken^ and several 
approache s 66 ' 67 indicate that there is a parity-breaking 
dimerization transition at a strong but finite value of 32, 
although it seems the answer may depend on the spe- 
cific lattice realization.— Thus, at least for some range of 
g = <?3 the coupling g 2 will flow to zero and one expects 
a phase transition for N = 3 on the way from mo = to 
mo = 00. For N > 4 the resulting CP N ~ 2 model with 
O = 7r remains dimerized, so no phase transition takes 
place. 
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FIG. 2: A sketch of the phase diagram of the model (|3j) 
in the vicinity of the AO-symmetric point 9 = —tt/2, 
in presence of two symmetry-breaking perturbations (|17p 
and (1431) : (a) the one-dimensional case with AT = 3; for 
N > 4 only the dimerized phase survives around the SU (AO- 
symmetrical point; (b-d) the two-dimensional case for dif- 
ferent values of the bare coupling g = nJ 1 Vl + A -1 ; here 
g^ — n 2 /[A(N — 1)]; for TV = 3 only, the degeneracy of the 
dimerized phase changes from fourfold at 9 < —tv/2 to twofold 
at 9 > -tt/2. 



m . The combined effect of the SU(N) h-> SU(N - 1)- 
breaking perturbation (|4*3")) and the >!?[/ (A^) i— > O(A^) one 
(TTT| is also transparent: taken together, those terms 
lower the symmetry to 0(N — 1), and for A/ > 4 the 
corresponding behavior as a function of 7 at finite mo 
can be inferred from the behavior of the model with 
N — > N — 1. In one dimension, for A/ = 3 and at 
large mo, 7 > favors a phase with dominant power- 
law XY-type nematic correlations (the XY2 phase in the 
classification of Schulz 7 -), while 7 < favors the Ising- 
type long-range antifcrromagnetic order. The transition 
from the XY-nematic to the dimerized phase is of the 
Berezinskii-Kosterlitz-Thouless type, and the transition 
from the AF-Ising to the dimerized phase belongs to the 
Ising universality class. The corresponding phase dia- 
grams are sketched in Fig. O 

For the spin-1 model {T]) that corresponds to N = 3, 
it is instructive to construct the effective Hamiltonian in 
the limit of strong single- ion anisotropy D> 1. Indeed, 
in that limit the Hilbert space at each site n is effectively 
reduced to the two spin-1 states |+), |— ), which can be 
identified with | |) and | J.) states of pseudospin-i. In 
the second order of perturbation theory in 1/D, the ef- 
fective Hamiltonian is given by the XXZ model in terms 
of pseudospin- i operators T n : 



L n,n' 
J xy 
J, 



n 

-J X y«r*, + rlrl ) + X«, , (46) 



-2 sin 6 



(cos 9 - sin eV2) 2 



2D 



J xy + 4cos#. 



One can see that for = —tt/2 the effective Hamiltonian 
is S'[/(2)-symmetric, in agreement with the continuum 
field description. Deviations from 9 = —tt/2 break this 
SU(2) symmetry, favoring AF or nematic order. 



SUMMARY 



In the two-dimensional usual poor-man's RG 

calculation yields the effective coupling g as a function of 
the anisotropy m§: 



!- ^[(JV-lJA-moarctan^]' 

Depending on the value of the bare coupling g = 
n^ 1 y/l + 1/A, there are three possible scenarios: (a) if 

g < <?jy = 7r 2 /[A(Y— 1)], then the system has long-range 
nematic or AF order all the way from mo = to mo = 00 ; 
(b) if gff < g < 9^_i, then the system is disordered (and 
dimerized) at mo = 0, but with increasing mo there is an 
ordering transition at mo ~ (27r/<7o)(l — 9/5^); finally, 

(c) 

if g > g N _i, the system stays dimerized at all values of 



We have studied the consequences of explicit symmetry 
breaking in the model of low-dimensional SU (N) antifcr- 
romagnet on a bipartite lattice, motivated by the physics 
of cold spinor bosonic atoms in optical lattices. Two pos- 
sible routes have been considered: lowering the SU(N) 
symmetry down to O(N) and to SU(N — 1). Physically, 
in cold atom systems those perturbations naturally arise 
due to the presence of the external magnetic field which 
controls the detuning from the Feshbach resonance and 
simultaneously causes the quadratic Zeeman effect. Both 
ways of the symmetry breaking result in rich sequences 
of transitions between dimerized, antiferromagnetic, and 
spin-nematic phases. The qualitative form of the phase 
diagram depending on the model parameters is estab- 
lished. It is shown that the physically interesting case 
N = 3 is special: perturbing the SU (3) symmetry leads 
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to nontrivial changes in the Berry phases, which are re- 
flected in the degeneracy of the dimerized phase. 

Acknowledgments 

I would like to thank Boris Ivanov, Ian McCulloch, 
Subir Sachdev and Kun Yang for stimulating discussions. 
The present work has been supported by the Heisen- 
berg Program Grant No. KO 2335/1-2 from Deutsche 
Forschungsgemeinschaft . 

APPENDIX A: EVEN-CHARGED SKYRMIONS 
AND MONOPOLES IN THE PERTURBED CP 2 
MODEL 

Consider a general q = 2 skyrmion solution of the un- 
perturbed (7 = 0) (l+l)-dimensional CP 2 model, which 
has the form 

Za = fa/\f\, fa = Cl a Z 2 + b a Z + C a . (Al) 

Now we would like to adjust the parameters of the above 
solution to satisfy z 2 = 0, making it suitable for 7 < 
0. Denoting the real and imaginary parts of the three- 
component complex vectors a, b, and casai, a 2 etc., we 
see that (ai, a 2 ) must be (up to a scale factor) a pair of 
mutually orthogonal unit vectors, and the same is true 
for (ci,c 2 ). We choose the coordinate system so that 
Ci || x and c 2 || y, and set b z = 2R > to fix the overall 
phase and norm. The following solution does the job: 

b = {C*,C,2R}, c = R 2 {l,i,0}, a = ai +ia 2 , 

ai = { - 1 - A 2 cos 2%, A 2 sin 2%, 2A sin x} , 

a 2 = { - A 2 sin 2%, 1 - A 2 cos 2%, -2A cos x} , (A2) 

where C is an arbitrary complex number and A, x are 
real. It is easy to convince oneself that this solution is 
nothing but the disguised Belavin-Polyakov skyrmion 79 
of the 0(3) nonlinear sigma model with the topological 



charge Q = 1. The correspondence between the CP 2 
field z and the sigma-model unit vector I is given by (. = 
—i(z* x z), and the 0(3) topological charge is determined 
by (j3"2"|) . One can easily see that the simplest Belavin- 
Polyakov solution (£1 + — £3) = Z/R translates 
into 

z= — {Z 2 + R 2 ,i(R 2 - Z 2 ),2iRZ\ 

V2(\Z\ 2 + R 2 ) X ^ h S 

which after a rotation Z 1— > Ze~ n l 2 becomes a special 
case of (HU) with A = 0, x = n/2, and C = 0. This so- 
lution describes a q — 2 skyrmion whose six constituents 
("zindons") sit at Z = ±R, Z = ±iR, Z = and Z = 00. 

In a similar way, one can show that in (2+1) dimen- 
sions a monopole of the CP 2 model with the even integer 
charge q = 2m, defined as a solution to the equation 80 

dA jL _qx }t 

£a ^ dx v ~ 2r 3 ' . { } 

where r 2 = a; 2 and the monopole is assumed to be at 

the origin, corresponds exactly to the hedgehog solution 

of the 0(3) model with a charge Q = m. Indeed, it is 

straightforward to check that the solution of the form 

j I cos 9 cos(m<p) — i sin(my) | 
z = —= < cos 9 sin(m<y9) -I- i cos(m</j) > , (A4) 
v2 - S - m Q J 

where 9 and tp are the spherical angular coordinates in 
the (2+l)-dimensional space, satisfies (|A3|) with q = 2m, 
satisfies z 2 = 0, and its corresponding 0(3) unit vector 
field £ = —i(z* x z) describes aQ = m hedgehog: 

I = { sin^ cos (p, sinfl simp, cos 0}. (A5) 

This confirms that the even-charged monopoles (|A4[) re- 
main exact solutions even in the perturbed case (but only 
for 7 < 0). Odd-charged monopoles are suppressed as ex- 
plained in Sect. IIII CI 
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